Topological effects at short antiferromagnetic Heisenberg chains 
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The manifestations of topological effects in finite antiferromagnetic Heisenberg chains is examined 
by density matrix renormalization group technique in this paper. We find that difference between 
integer and half-integer spin chains shows up in ground state energy per site when length of spin 
chain is longer than ~ £, where £ ~ exp(7rS) is a spin-spin correlation length, for spin magnitude 
S up to 5/2. For open chains with spin magnitudes S = 5/2to5' = 5, we verify that end states 
with fractional spin quantum numbers S' exist and are visible even when the chain length is much 
smaller than the correlation length f. The end states manifest themselves in the structure of the 
low energy excitation spectrum. 



PACS Numbers: 75.fO.Jm, 75.40.Mg 

The importance of topological term in antiferromag- 
netic spin chains was first pointed out by Haldane by 
mapping-|the Heisenberg model onto nonlinear er-model 
(NLcrM)BH. As a direct manifestation of the topological 
effect, Haldane pointed out that (infinite) integer spin 
chains have a gap in their excitation spectrum, whereas 
half-integer spin chains are gapless. The gap A in integer 
spin chains scale as A ~ Je~ 7rS , where J is the exchange 
coupling, and S the spin magnitude. Topological effect is 
expected to manifest itself also at properties not directly 
related to excitation spectrum, for example, oscillation in 



state energy per site for antiferromagnetic Heisenberg 
open chains with Hamiltonian 
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L-l 
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grouncLstate energies as ' 
integercl. More recently Ng 



..changes from half-integer ta 
I pointed out that end statescl 



appear as lowest energy excitations in long open spin 
chains due to the same topological term, and the struc- 



where the chain length is L and S is spin operator. First 
we consider the large L limit. In this limit the ground en- 
ergy 60(6*) can be obtained accurately as half of the differ- 
ence of the energies for length L and L+2 chain in DMRG 
methodd. We keep m = 500 states in DMRG, and the 
biggest truncation error is 10 -6 . We obtain the following 
converging results for thermodynamic limit from data of 
L up to several thousands: 



ture of end states spectrum is_ determined by the spin 
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magnitude S of the spin chainO. In this paper we shall 










study numerically the manifestation of topological effect 
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rable with the lowest spinwave excitation energy of the 
finite spin chain and becomes unobservable. If this is the 
case, then we expect that the topological effect will be 
hard to observe for spin chains with large value of spin 
magnitude S, where the correlation lengths are very long. 
However, the derivation of the NLerM assumes only that 
the length scale is much larger than lattice spacing and 
the topological term exists as long L » 1, independent 
of the correlation length £. Therefore, we expect that 
topological effect may manifest itself in some properties 
of spin chain even when chain length L is less than £, as 
long as L » 1. In this paper, we shall study two prop- 
erties of finite spin chains: the ground state energy and 
the end states. We find that oscillation in ground state 
energies between integer and half-integer spin chains is 
indeed observed numerically when chain length L > £. 
However the topological end states are much more ro- 
bust and appear in finite spin chains even when L << £. 
We use the DMRG methods to calculate the ground 



The ground state energy for S = 1/2 chain is exactly 

known, with e (0.5) = 1/4 - ln(2) 0.4431471. The 

energy for S q 1 chain has also been obtained to very 
high accuracy^, where e (l) = -1.401484038971. In a 
paper demonstrating the k =_1 SU(2) WZW low en- 
ergy behavior of S — 3/2 chain3, the ground state for it 
has been accurately obtained, with eo(1.5) = —2.82833. 
The energy for S — 2 chain was obtained in another 
paper demonstrating its massive relativistic low energy 
propertyQ, eo(2) = —4.761244. Our calculation is in 
agreement with these earlier studiesEJ. 

The energies can be compared with a 1/S expansion: 



e (5) = S 2 + (2/tt - l)S + a Q + aJS + . . 



(2) 



where ±he first two terms were obtained from spin- wave 
theoryB In Fig. | we plot -S 2 - (1 - 2/tt)5 - e (S) 
as a function of 1/S. Our numerical result shows clearly 
the oscillatory nature of eo(5) between integer and half 
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integer spin chains and eo(S) cannot be fitted by a single 
mono-tonic function of 1 / S. Assuming that the oscillatory 
behaviour is coming from topological terms in NLaM 
we expect that the oscillatory part of the ground state 
energy scales as exp(-irS) and is nonanalytic in a l/S 
expansion. This is indeed consistent with our numerical 
result as can be seen from the fitting in Fig. [|. 
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Finite size scaling for ground state 
for antiferromagnetic Heisenberg chains. 
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Fig.1 Lou et al, PRL 
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FIG. 1. Ground state energy for antiferromagnetic Heisen- 
berg chains. —S 2 — (1 — 2/-k)S — eo(S) vis l/S is plotted and 
fit by polynomial of l/S. The ground state energies for inte- 
ger spin chains and half integer spin chains cannot be fitted 
by a single monotonic curve. The difference decreases to zero 
very fast as S increases. We note that the 5=1/2 point is 
not used in the fitting process 



Fig. 2 Lou et al, PRL 




(eo — EO/L)/^/ S(S + 1)) versus inverse of the renormalized 
length L/£ is plotted for 5 = 1/2, 1, 3/2, 2, and 5/2 with ap- 
proximated values of £ = 2, 6.3, 19, 49, and 140, respectively. 
Integer S and half-integer S chains behave differently in fi- 
nite size scaling, but the difference vanishes for small L when 

i/L > i, 

As is shown in Fig. [t], for S > 3, the energy differences 
between integer and half integer spin chains due to topo- 
logical effect are very weak. To examine the chain length 
dependence of ground state energy we plot in Fig. ^ the 
the length dependence of the finite length correctio n of 
the ground state energy L\e§L — E0(L)]/ ^/S(S + 1) ver- 
sus the inverse of renormalized length (i/^) _1 for various 
value of spin magnitude S up to S = 5/2. The correla- 
tion lengths are estimated independently through various 
fittings and the scaling behaviour is insensitive to fluc- 
tuations in values of £ within ~ 10 percent. EQ(L) is 
computed using DMRG method for chains with periodic 
boundary condition (PBC). It is clear from the figure 
that the ground state energies of integer and half-integer 
spin chains behave differently in finite size scaling, and 
the difference vanishes as L < £, in agreement with naive 
expectation based on Haldane &a/p argument. 

It was pointed out in Ref. p] that topological effect 
also leads to appearance of end states in open spin chains. 
For integer spin chains with length L > > £ end spins with 
magnitude S/2 appear at each end of spin chain and are 
coupled to each other with effective Hamiltonian 



Heff — J(L)S i • S'l, 



(3) 



where S'i and S'l are the two end spins on the ends of 
the open chain. J(L) ~ Je^ L ^ is positive for even L 
and negative for odd L. For half-integer spin chains end 
spins with magnitude (5 — 1/2) /2 appear and are coupled 
to each other with same effective Hamiltonian (3) except 
that J(L) ~ L\'l(L) f° r l ar S e -k- Notice that bulk spin- 
wave excitations have energies scale as ~ 1/L for large L 
and the lowest spin excitations of both integer and half- 
integer open spin chains are end-spin excitationSjidicn L 
is large enough. Note that end states for S = 1LIE3 and 
S = 3/2, 2lLj open chains have been observed numeri- 
cally in long chains L is > £ with properties in agreement 
with end state theoryQ. The question of interests here 
is whether these end states remain robust as spin value 
S increases, with length of spin chains reduce to L < £. 
We expect that the end states may stay robust because in 
general energy level crossings have to occur if end states 
are moved out of the low energy excitation spectrum. 
We shall show in the following that end states are ob- 
served as lowest energy excitations in open spin chains 
for large spin magnitude S up to 5 with chain lengths 
much smaller than £. 

We consider spin chains with spin magnitudes 5/2 < 
S < 5 and with even number of sites with chain lengths 
from L = 4 to 30. Notice that by restricting ourselves 



2 



to chains with large value of 5 and short lengths < 30 
we are always in the limit L < £ in our numerical study. 
By keeping m = 450 states in DMRG. We calculate the 
lowest energy states in sectors with 5* ot = 0,1,. ..,5+1. 
The biggest truncation error is 10~ 4 for 5 = 5. Accord- 
ing to end state theory, the lowest 25' + 1 states in the 



excitation spectrum should have 5 = 0, 1, 



25' 



increasing order of energy, where 5' = 5/2 for integer 
spin chains and S' = (S — l/2)/2 for half-integer chains, 
the bulk spinwave spectrum appears only above these 
states. 



Fig.3 Lou et al, PRL, integer S 
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FIG. 3. Low energy excitations times chain length for 
integer spin 5 = 3 to 5 open chains. The ground state is of 
total spin S tot = 0. The 5 + 1 lowest excited states from low 
energy to high energy are S tot = 1, . . . , 5, and S tot = 1 from 
bottom to top. 



Fig. 4 Lou et al, PRL, half-integer S 
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We show the low energy excitation spectrum for open 
integer spin chains in Fig. 3. Starting from the ground 
state energy EO, we label the i th excited states' energy 
as Ei. For integer spin chains, we find indeed that the 
lowest 5+1 energy levels are ordered in sequence with 
gtot _ o, 1, . . . , 5 from ground state to the S th excited 
states. Above this series of states, spinwave excitations 
appear, starting with total spin S tot = 1. The corre- 
sponding low energy spectrum for open half integer spin 
chains is shown in Fig. 4, and the same structure is 
observed- the lowest 5+1/2 energy levels are ordered in 
sequence with S tot = 0, 1, . . . , 5 — 1/2 from ground state 
to the (5 — l/2) th excited states. Above these series of 
states, spinwave excitations appear. These are exactly 
the spectra predicted by the end state theory. 

To examine the nature of these low lying states further 
we investigate the validity of Eq. (S) in describing the en- 
ergies of these states. It is easy to show that Eq. (|J) pre- 
dicts that the number series y(i) — (Ei — EQ) /(El — EO) 
for the end states is given by the simple mathematical ex- 
pression y(i) = i(i + l)/2, independent of chain lengths 
L and spin magnitude 5. 
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Fig. 5 Lou et al, PRL S=3,4,5 
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FIG. 5. Edge excitation energy for integer spin chains. 
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FIG. 4. Low energy excitations times chain length for 
half-integer spin S = 5/2 to S = 9/2 open chains. The 
ground state is of total spin S tot = 0. The 5-1/2 
lowest excited states from low energy to high energy are 
5' ot = l,...,5-l/2. 
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Fig. 6 Lou et al, PRL S=5/2,7/2,9/2 
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FIG. 6. Edge excitation energy for half integer spin chains. 

In Fig. 5 and Fig. 6, we plot respectively, y(i) for the 
25" lowest energy levels given by DMRG versus 1/ ln L 
for integer spin chain and half integer spin chains. It is 
clear that Eq. (||) describes quite well the qualitative be- 
haviour of the lowest 2S" states of the excitation spectra 
for all chains with different spin values and lengths under 
investigation, despite that we are already in the L < £ 
limit. 

Summarizing, we study in this paper using DMRG 
method several manifestions of topological effect in fi- 
nite Heisenberg spin chains. We confirm the oscillation 
of ground state energies between integer and half-integer 
spin chains. The magnitude of the oscillatory part of 
ground state energies goes down very rapidly as S in- 
creases and is consistent with an e~ vS behaviour. The 
oscillation becomes unobservable when chain length L de- 
creases to less than £, in agreement with argument based 
on Haldanc gap. We have also study open spin chains 
where topological effect also manifest itself as end states. 
Surprisingly, we find strong numerical evidence that end 
states stay robust as chain length L's decrease and are 
observable even when L << £. Our result confirms the 
theoretical expectation that the topological effect is a ro- 
bust property of quantum spin chains and exists as long 
as chain length L » 1, independent of the appearance 
of Haldane gap in the excitation spectrum. 
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